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Abstract. Let Q be a bounded smooth domain in R N . We con- 
sider the problem u t = Au + V(x)u p in fi x [0, T), with Dirich- 
let boundary conditions u = on dfl x [0, T) and initial datum 
u(x, 0) = Mip(x) where M > 0, tp is positive and compatible with 
the boundary condition. We give estimates for the blow up time 
of solutions for large values of M. As a consequence of these esti- 
mates we find that, for M large, the blow up set concentrates near 
the points where ip p ~ l V attains its maximum. 



1. INTRODUCTION 

In this paper we study the blow-up phenomena for the following 
semilinear parabolic problem with a potential 

u t = Au + V{x) u p inQx(0,T), 

(1.1) u(x,t) = on dfl x (0,T), 
u(x, 0) = Mip(x) in a 

First, let us state our basic assumptions. They are: Q is a bounded, 
convex, smooth domain in M> N and the exponent p is subcritical, that 
is, 1 < p < (N + 2)/(N — 2). The potential V is Lipschitz continuous 
and there exists a constant c > such that V(x) > c for all x G f2. As 
for the initial condition we assume that M > and that (p is a smooth 
positive function compatible with the boundary condition. Moreover, 
we impose that 

min V{x) 

(1.2) MAip + ^ MV > 0. 
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We note that (jl.2j) holds for M large if Ap is nonnegative in a neigh- 
borhood of the set where <p vanishes. 

It is known that, and we will prove it later for the sake of complete- 
ness, once ip is fixed the solution to (jl.lj) blows up in finite time for any 
M sufficiently large. By this we understand that there exists a time 
T = T(M) such that u is defined in Q x [0, T) and 

lim \\u(-,t) \\ L oo (n) = +00. 

The study of the blow-up phenomena for parabolic equations and 
systems has attracted considerable attention in recent years, see for 
example, jB], jBB], pO] . [HK2] . [HVj . jHVT] . jHV2] . [M], [Z] and the 
corresponding references. A good review in the topic can be found in 
|GV2j . When a large or small diffusion is considered, see |lYj . |MYj . 

Important issues in a blow-up problem are to obtain estimates for 
the blow-up time, T(M), and determine the spatial structure of the set 
where the solution becomes unbounded, that is, the blow-up set. More 
precisely, the blow-up set of a solution u that blows up at time T is 
defined as 

B{u) = {x/ there exist x n — > x, t n /* T, with u(x n ,t n ) — > 00}. 

The problem of estimating the blow-up time and the description and 
location of the blow-up set has proved to be a subtle problem and has 
been addressed by several authors. See for example [SGKMj, |GV2j 
and the corresponding bibliographies. 

Our interest here is the description of the asymptotic behavior of 
the blow-up time, T(M), and of the blow-up set, B(u), as M — > 00. It 
turns out that their asymptotics depend on a combination of the shape 
of both the initial condition, (p, and the potential V. Roughly speaking 
one expects that if (p = 1 then the blow-up set should concentrate near 
the points where V attains its maximum. On the other hand if V = 1 
the blow-up set should be near the points where p attains its maximum. 
Just to see what to expect, if we drop the laplacian, we get the ODE 
u t = V(x)u p with initial condition u(x, 0) = M<p(x). Here x plays the 
role of a parameter. Direct integration gives u(x, t) = C(T — i)^ 1 /^ 1 ) 
with 

M 1_p 

T = . 

(p- l)V(x)ipP~ 1 (x) 

Hence, blow-up takes place at points xq that satisfy V(xo)p p ~ 1 (xo) = 
max x V(x)p p ~ 1 . Therefore, we expect that the quantity that plays a 
major role is (max x V(x)p p ^ 1 (x)). 
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Theorem 1.1. There exists M > such that if M > M the solution 
°f (ll-lj) blows up in a finite time that we denote by T(M). Moreover, 
let 

A = A(ip, V) := j— 

then there exist two positive constants C\, C 2 , such that, for M large 
enough, 



;i.3) --^r < T(M)M^ --^r< -^r, 



and the blow-up set verifies, 

(1.4) ^-\a)V{a) > i - /or a// a G 

where 7 = min{^-, |). 

Note that this result implies that 



lim T(M)M 



v-i 



A 



M-+oc p — 1 

Moreover, it provides precise lower and upper bounds on the difference 
T(M)M p - 1 - 

We also observe that (jl.4j) shows that the set of blow-up points con- 
centrates for large M near the set where ip p ~ 1 V attains its maximum. 

If in addition the potential V and the initial datum ip are such that 
tp p ~ 1 V has a unique non degenerate maximum at a point a, then there 
exist constants c > and d > such that 

^(a^a) - if- l {x)V{x) > c\d - x| 2 for all x e B(a, d). 

Therefore, according to our result, if M is large enough one has 

C 

la — a\ < — ^ for any a G B(u), 

with 7 = min( 2 ^-, |). 

Throughout the paper we will denote by C a constant that does not 
depends on the relevant parameters involved but may change at each 
step. 

2. Proof of Theorem fTTTl 

We begin with a lemma that provides us with an upper estimate 
of the blow-up time. This upper estimate gives the upper bound for 
T(M)M p ~ l in (|1.3J) and will be crucial in the rest of the proof of The- 
orem 
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Lemma 2.1. There exist a constant C > and Mq > such that for 
every M > M , the solution of (jl.lj) blows up in a finite time that 
verifies 

(2.1) T(M)< £ - + ^ . 

Proof: Let a G O be such that 

tf-\a)V(a) = ma^ V p -\x)V(x), 

X 

L the constant of Lipschitz continuity of V, and K an upper bound for 
the first derivatives of if and L. 

In order to get the upper estimate let M be fixed and e = e(M) > 
to be defined latter, small enough so all functions involved are well 
defined. Pick 

6 = ^ 

then 

V(x) > V{d) — - and > f(a) — e for all x G B(a, 5). 

Let w be the solution of 

w t = Aw+ (V(a) - |) w p in B(a, 5) x (0, T w ), 

w = on dB(a,S) x (0,T W ), 

w(x, 0) =M((p(a) -e), in B(a,S) 

and T w its corresponding blow up time. A comparison argument shows 
that u > w in B(d, 5) x (0, T) and hence 

T<T W . 

Our task now is to estimate T w for large values of M. To this end, 
let Ai(<5) be the first eigenvalue of —A in B(d,S) and let <pi be the 
corresponding positive eigenfunction normalized so that 

4>i(x) dx = 1. 

B(a,S) 



That is, 



Now, set 



= Ai(<J)0i, mB(a,8), 
'i = on dB(d, 5). 



= / w(x,t)(j) 1 (x)dx. 

'B(a,S) 
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Then $(t) satisfies $(0) = M(cp(a) - e) and 



$'(£) = / w t (x,t)(j) 1 (x)dx 

'.8(5,(5) 



Au7(x,t)0i(a;) + - |) w p (a;, t)0i(a?) ) dx 



8(5,(5) 

> — Ai(5) / w{x,t)<pi{x) dx 

JB(a,S) 

+ (v(a) — -j ( / w(x,t)(f)i(x) dx 

^ 2/ \JB(a,8) / 

= -AxOWf) + (v(a) - 1) 

Let us recall that there exists a constant D, depending on the di- 
mension only, such that the eigenvalues of the laplacian scale according 
to the rule Ai(5) = D5' 2 . 

Now, we choose e such that 

Ai(tf) = D5- 2 = D = ~(M(<p(d) - e))^ 1 . 

So, e is of order 

C 



p-i 



Choose Mo such that for M > Mq the resulting e is small enough. 
Then for any M > Mo we have that 

(2.2) $'(t)>(V(a) -e)$(t)* 

for all t > for which $ is defined. 

Since $(0) = M(ip(a) — e) and T w is less or equal than the blow up 
time of $ integrating (|2.2|) it follows that 

T < I 

w ~ MP' 1 (p - 1) (V (a) - e) (p(a) - e)?- 1 

< ' ... .+ c 



Mp- 1 ^ - 1)^(0)^(0)^ ' M^Mf" 1 ' 
for all M > M . □ 

Now we prove a lemma that provides us with an upper bound for 
the blow up rate. We observe that this is the only place where we use 
hypothesis ()1.2j) . 

Lemma 2.2. Assume (|1.2jl . Then there exists a constant C indepen- 
dent of M such that 

u(x,t) < C{T-t)~p^. 
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Proof: Let m = minV. Following ideas of ll'Mcj. set 



m 

V = Ut — —u- 



Then v verifies 

v t - Av - V{x)pu p ' 1 v = —p(p - l)u p ~ 2 \ Vu\ 2 > in Q x (0, T), 
v = onfflx(0,T), 

v(x, 0) = MA(p + (V(x) - M p ipP > in a 

Therefore v > and hence 

^ m p 
2 

Integrating this inequality from to T we get 

u(x,t) < — = C(T -ty~\ 

(m(p- l)(T-t)))~ 

as we wanted to prove. □ 

We are now in a position to prove Theorem 11.11 

Proof of Theorem ll.lt The idea of the proof is to combine the 
estimate of the blow-up time proved in Lemma 12.11 with local energy 
estimates near a blow-up point a, like the ones considered in |GKlj and 
[GK2 , to obtain an inequality that forces ip p ~ 1 {a)V{a) to be close to 
maXj ip p ~ 1 V. 

Let us now proceed with the proof of the estimates on the blow-up 
set. We fix for the moment M large enough such that u blows up in 
finite time T = T(M) and let a = a(M) be a blow up point. As in 
GK2], for this fixed a we define 



w{y,s) = (T - t)p-^u{a + y(T - t)a ,t)\ t=T ^ e -sy 
Then w satisfies 

11 _s 

(2.3) w s = Aw y-Vw w + V(a + yTe 2 )w p , 

2 p 1 

in U se (o,oo)^( s ) x {s} where f2(s) = fi a (s) = {y : a + yTe~i E f2} with 
w(y, 0) = Tv^ip^a + yT^). The above equation can rewritten as 

11 _s 

w s = - V(pVw) w + Via + yTe 2 )w p 

P P - 1 

where p(y) = exp(^j^). 
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Consider the energy associated with the "frozen" potential 

V = V(a), 

that is 

E(w)= [ (\\Vw\ 2 + -^w 2 - ^V{a)w?A p{y)dy. 
Ms) \2 2(p- 1) p + 1 / 

Then, using the fact that is convex, we get 

~T<- {ws?p{y) dy + / (V(a + yTe-i) - V{a))w p w s p{y) dy. 
as Jn(s) Jn(s) 

Since V(x) is Lipschitz and w is bounded due to Lemma 12.21 then 
there exists a constant C depending only on N, p and V, recall that 
the constant in Lemma [2.21 does not depend on M, such that 

1/2 

[w s ) 2 p(y) dy 



dE 
ds 



<-J(w s ) 2 p(y) dy + Ce-iT(J( 



Maximizing the right hand side of the above expression with respect 
to f (w s ) 2 p(y) dy we obtain 

— < Ce~ s T 2 
ds 

and integrating is s we get 

(2.4) E(w) < E(w ) + CT 2 . 



Since w is bounded and satisfies (|2.3|) . following the arguments given 
in |GKlj and |GK2j . one can prove that w converges as s —>■ oo to a 
non trivial bounded stationary solution of the limit equation 

(2.5) = Az--yVz —z + V(a)z p 

2 p — 1 

in the whole Mr. 

Again by the results of |GKlj and |GK2j . since p is subcritical, 1 < 
p < (N + 2)/(N — 2), the only non trivial bounded positive solution of 

(12. 5j) with V(a) = 1 is the constant (p — l)^^. A scaling argument 
gives that the only non trivial bounded positive solution of ()2.5|) is the 

constant k = k(a) given by 

k(a) = — . 

(V(a)(p-l))—i 

Therefore, we conclude that 

lim w = k(a) 
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if a is a blow-up point. Also by the results of |GKlj . |GK2| we have 

(2.6) E(w(; s)) -> E(k(a)) as s -> oo, 

where 

E(k(a)) = J f-J—^k(a)f - J—V(a)(k(a)rA p(y)dy 

By (12. 4j) and ()2.6|) we obtain that, if a is a blow-up point, then 
E(k(a)) < E(w ) + CT 2 . 

where w (y) = w(y, 0) = Tp 31 Mip{a + yT^). 

As if is smooth, yp(y) integrable, and Tp^M is bounded by Lemma 
12.11 there are constants C independent of a such that for M > Mo 

E(w(-,0)) = [ (hvw (y)\ 2 + 



n(o) 



2 



>(p- !) w q^) J p(y) d y 



I [^-V{a)wl + \y))p{y)dy 
'o(o) VP + 1 / 

< / fi(T^iM) 2 T|V^(a)| 2 ')p( 2 /)^ 
Jn(o) \ z J 

+ LS^) {T ^ Mv{a)) ^ p{y) ^ 

- f (^V{a){T^Mtp{a)yA p(y) dy 

Mo) \p + 1 / 
+crl + cri 

Therefore, since |V</?| is bounded, 

E(w(-,0)) < [ (^^(T^M<p(a))Ap(y)dy 

' V(a)(T^M^(a)) p+1 )p(y)dy 



O(0) 



'n(o) VP + 1 
+CT§ + CT§. 

Or, since T < 1 for M large 

0)) < E(T^M<p(a)) + CT^. 

Hence we arrive to the following bound for E(k(a)) 

(2.7) E(k(a)) < E(w(-, 0)) + CT 2 < E(T**M<p(a)) + CT^ . 
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Observe that if b is a constant then the energy can be written as 

E(b) = VF(b), 

where T is the constant 

r = J p(y) dy 

and F is the function 

F(z) = ( , l , z 2 —V(a)z p+1 

1 ; \2{p-l) p+1 V ; 

As F attains a unique maximum at k(a) and F"(k(a)) = — 1 there 
are a and (3 such that if \z — k(a)\ < a then 

*"M < ~, 

and if \F(z) - F(k(a))\ < (3 then 

\z — k(a)\ < a. 

From (J2.7|) we obtain 

F(k(a)) < F(T^M<p(a)) + CT^ . 
If Mi is such that C{T{M 1 ))^ = (3 then for M > max(M , Mi) 

/3 > CT5 > F(Jfc(a)) - F(T^M(p(a)). 
Hence by the properties of F, 

\k(a) -T^M<p(a)\ < a. 

Therefore 

CT^ > F{k{a)) - F{T^M<p{a)) > ^{T^M<p{a) - jfc(a)) 2 . 
So, using Lemma EU 

k(a) - CT\ < T^M<p(a) 

. <fip) , Cip{a) 

(2.8) —j- i — — + 



(p _ i)Fiyi=i (a>(a) Ms 
k{a)t){a) 



Ml ' 



where 



0(a) 
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and a is such that 

ip p -\a)V(a) = max<f?-\x)V(x). 

X 

Recall that 



Mp- 1 ' 
Therefore, we get 

k(a)(l-6(a))<-^ T l + -— T < 



with 7 = min( p -j^- - 



4 ' 

As V is bounded we have that k(a) is bounded from below, hence 

(l-6(a))< ( 



that is, 

C 



0(a) > 1 



and we finally obtain 

(2.9) <p(a)V(a)** > <p(a)V(a)^ 



Mi 

This proves (jl.4j) . 

To obtain the lower estimate for the blow-up time observe that from 
(I2~T1 and the fact that V(a) > c> we get 

Vp(a)y(a)^ C 

p(a) > 



V(a)~ V(a)^Mf 



(2.10) ^ ^(a)F(o)5=i c 



> C > 0. 



Inequality ()2.8|) gives us 



1 — - CT3 < T^M V {a) 



(V(a)(p-l)y- 
Hence 

1 CT? 



< T~M. 



y{a){V{a){p-l))^ <^( a ) 
By fl2~TUJ) and ifp-\a)V(a) = max, ^' l (x)V{x) we get 

1 i _ i . 
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and using 



~ MP- 1 
we obtain 



C 

< T^M 



V {a){V{a)(p-l))— M— 
as we wanted to prove. □ 
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